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Abstract 

We study the time evolution of the radial perturbation for self-gravitating soliton and black- 
^^ \ hole solutions in a generalized Skyrme model in which a dilaton is present. The background 

solutions were obtained recently by some of the authors. For both the solitons and the black 
holes two branches of solutions exist which merge at some critical value of the corresponding 
parameter. The results show that, similar to the case without a scalar field, one of the branches 
is stable against radial perturbations and the other is unstable. The conclusions for the linear 
stability of the black holes in the generalized Skyrme model are also in agreement with the results 
from the thermodynamical stability analysis based on the turning point method. 



1 Introduction 

Our intuition for the properties of the solutions describing self-gravitating objects in general rela- 
tivity is based, to a large extent, on some exact solutions which belong to the Kerr-Newman class 
C^ \ of black holes. For these solutions uniqueness theorems, and theorems stating that globally regular 

self-gravitating solutions (solitons) do not exist, have been proven rigorously in the case of vacuum or 
linear matter models such as Maxwell electrodynamics [I]-[5]. As the investigations in the last two 
decades revealed, the standard intuition often fails when nonlinear matter models are considered [B], 
which makes the study of self-gravitating solutions in such models vital for fundamental physics. 

One of the effective nonlinear matter models which has attracted much attention is Skyrme's 
theory [7], [5]. In this theory baryons are described as solitons in an effective theory of mesons. The 
interest in Skyrme theory was revived in the 1980s when it was found that the Skyrme Lagrangian 
can be derived from quantum chromodynamics (QCD) in the low-energy regime. 

Self-gravitating solutions in Skyrme theory were considered for the first time by Luckock et al. |9| . 
The solutions in Einstein-Skyrme (ES) theory are nonunique, and those with a nontrivial Skyrme field 
can be divided into two branches. The first branch of solutions has a well-defined flat-space limit. It 
was obtained by Droz, Heusler and Straumann [10]. The authors found that these solutions are stable 
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against spherically symmetric perturbations [11] [12] . The second branch of solutions was discovered 
by Bizon and Chmaj soon after that J13j . This branch has no flat-space limit and it is unstable, as the 
authors' analysis revealed. The stability of the ES solitons has also been studied in [HI [151 [16] . There 
is also a branch of solutions that has a trivial Skyrme field and coincides with the pure Schwarzschild 
black hole. Self-gravitating solutions in Skyrme theory, both black holes and solitons, have also been 
studied in a series of papers [T7]-[27]. 

Different modifications of Skyrme theory have been considered in order to cure some of its de- 
ficiencies which are present in the original version of the theory 7 . One possible generalization is 
the inclusion of a dilaton. The dilaton is added in the theory to restore scale invariance which is 
also characteristic for the underlying QCD. It has also been considered as a source of additional 
intermediate-range attractive forces which are vital for the formation of stable multisoliton configura- 
tions such as nuclei and baryon stars. A Generalized Skyrme Model (GSM) which includes a dilaton 
has been derived from QCD in the low-energy regime in [28] [29l [3^ . In a recent paper "ST we re- 
ported numerical solutions describing self-gravitating solitons and black holes in the GSM. They are 
generalizations of the soliton and black-hole solutions that have been obtained numerically in |10| and 

m- 

The aim of the current paper is to study the response of the self-gravitating GSM solutions [3T] . 
both soliton and black-hole types, to small radial perturbations and, in particular, to determine if the 
inclusion of the dilaton in Skyrme theory changes the stability properties. We study the quasinormal 
modes (QNMs) of the solutions by evolving the time-dependent wave equations. 

The problem of studying the QNMs and the stability of the GSM solutions is mathematically 
more complex than that of the ES solutions since in the former case a system of two coupled wave 
equations for the perturbations of the Skyrme field and the dilaton has to be solved even though the 
considerations are restricted to radial perturbations (in the ES case the problem is reduced to only 
one wave equation for the Skyrme field). What makes the problem even more difficult is that the wave 
equation for the perturbations of the dilaton contains a potential which is not vanishing at infinity, 
i.e. the scalar field is massive and the time evolution of the perturbations has some specific properties 
[32]-[38]. 

The paper is organized as follows. The GSM coupled to gravity is briefly presented in Section [2] 
In this section the time-dependent field equations are given. The system of coupled equations for the 
radial part of the perturbations of the Skyrme held and the dilaton is derived in Section [3] and solved 
with the proper boundary conditions numerically in Section [J] In Section [5] a summary of the results 
is given. 

2 The Generalized Skyrme Model 

Action 

Let us briefly introduce the model considered in [3T]. We start with the following action: 

5./<i..V-j(-^ + L„). (1) 

The flat-space Lagrangian of the GSM can be found in [30]. When gravity is included the GSM 
Lagrangian is naturally generalized to the form 

Lm = ^/' exp(-2a) Tt[V^UV^U+] + ^^ exp(-2a) g^'^'d^ad^a (2) 

^ Tr[(V^C/)C/+,(V.L/)[/+]2 + yGSM(fT), 



32e2 



where the derivatives have been substituted with covariant derivatives. Here U is the SU(2) chiral 
field, a is the dilaton, V^ is the covariant derivative with respect to the metric g^^, /^r is the pion 
decay constant, e is the Skyrme constant, Cg is the gluon condensate, Nf is the number of flavors, and 
e = 8A^//(33 — 2Nf). The first two terms in ([2]) are the kinetic terms for the chiral and the dilaton 



fields. The third term is the one introduced by Skyrme for the stabilization of the soliton solutions. 
The potential of the dilaton field is given by 



Vgsm(o') = 



48 



exp(— 4ct) — 1 H — (1 — exp(— £cr)) 



(3) 



The dilaton couples only to those terms of Lagrangian density that break the scale invariancqj. 
Instead of a it is more convenient to work with the function $ which is defined by 

$ — exp(— cr). 



(4) 



Reduced Lagrangian 

We are going to restrict our considerations to the spherically symmetric case. In [ST the hedgehog 
ansatz for the chiral field 

U = exp[r • YF[r, i)] (5) 

was chosen. Here r are the Pauli matrices and f is a unit radial vector. With the following time- 
dependent ansatz for the metric, 



the Lagrangian ^ takes the forn|j 
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a 



where 



u = a;^$2 + 2 sin^ F, v=[ 2$^ 



~.^. IQnGh^^ .. -fNb 
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sin^ F, 



$ 



e 



We have introduced the following constants: 



a^SnGfl b^SnG^, 7 = 2^G^, N=^, 



(6) 

(7) 

(8) 
(9) 

(10) 



and dimensionless variables r — ef^^t, x — efj^r. The derivative with respect to the dimensionless 
time coordinate r is denoted by a dot, while the derivative with respect to the dimensionless radial 
coordinate a; is denoted by a prime. Below we will also use the parametei|2| 



Drff = 



jN 



2 f 2 ■ 



2 a e2/2 
For the number of flavors, we fixed the value Nf = 2, so A^ = 1 



(11) 



Time-dependent field equations 

The Einstein equations have the following form: 

Gni, = 



-It 



(12) 



^For more details we refer the reader to |31) . 

^The notation we choose here is sUghtly different from that in '31'. It faciUtates the comparison of the equations 
and the results to the ES case 1101 — 1131 . 

•^The parameter a is two times bigger that the parameter a used in |10 | — 113| and the parameter -Dcff is chosen to be 
the same as in | 30| . 
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The (tt), (rr), and (tr) components of (fT^ are 



[e " (1 - xa') - l] = -a ( -uw + -v + -x'^z ) + -x'^V, 
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The combination of equations ([H|) and p^ gives the following useful expression: 

I 1 - -av + -x^V I - 
X " " 



— — 1 av H — X V 
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The time-dependent field equations for F and <f> obtained from ([7]) are 



a-X 



ui^ + {uF) 



x'-«' 



uF' + {uF' 



(13) 

(14) 

(15) 
(16) 

(17) 
(18) 

(19) 



^UF {e"-^F^ -F'^)- ie%f = 0, (20) 



a<^-X 



^-^x^^ + {x^^) 



^^^x2$' + (x2$')' 



+ ^u$ (e"-^F^ -F'^)- ^e"w$ + ^^e"14 = 0, (21) 
27V V J 2N 2aN 

where (..)$ denotes the partial derivative with respect to $, and {■.)f denotes the partial derivative 
with respect to F. 



3 Equations for the radial perturbations 

We reduce our considerations to radial perturbations 

a(T, x) — ao{x) + da{T, x), 

x{t,x) ^Xo{x) + Sx{t,x), 

FiT,x)^Fa{x)+dF{T,x), 

$(r, x) = $o(a;) + S^{t, x), 

and follow the scheme presented in [TT] . It turns out that the evolution of the Skyrme field and the 
scalar field perturbations, 6F and 5^, respectively, can be studied independently from the perturba- 
tions of the metric. The equation for 6F, obtained from (P0|) . is 



^""-''''uqSF^uoSF" • ' -^^ "° 



uo + Uo]SF' + uq^Fq (5$ '+ 



Un*Fn ,5$ + ( Fr': + ^° ^"° F^j Su + u'^pF^5F - ]^F^''5uF- 



^0$^ ' 







ie"«,S«j. - ie"°t;oF^a + »o J^o "^^ ' /" ' ■ (22) 



From (pTj) we obtain the following equation for the perturbations of the scalar field S^ 



I 



V 2 /TV 

\ , 1 / a;2 ^ 

2iV 2A^ V a 

1 / x^ ~ \ b\' - ba' 
-^e"« -i;o* + —1^0* fc + 2:'$o-^^^ , (23) 

where 

bu = uqfSF + uo$(5$, bv = vqfSF + vo^b^, 

bup — UQppbF + u^p<i,b<^, bvp — voppbF + VQp^b^, 

bu^ = uo$_F'5F + uo$$(5$, bv,p = v^^pbF + i;o$$(^<f>, (5V$ = Vo$$(5$. 

Throughout the paper, the lower index (..)o means that the corresponding quantity refers to the 
background static solution. Lower indices F and $ denote the corresponding partial derivatives. The 
variation of eq. (J16p gives 

ba=- (uoFobF + x^N<^Qb'^] . (24) 

The integration of the above expression with respect to r gives 

ba^- (uoFobF + x^M'ob^) (25) 

and it allows us to relate the perturbations of the metric functions ba to the perturbations of the 
matter fields bF and (5$. Another useful relation can be obtained from ([T5)) 



bx' -Sa' _ e"o 
2 ~ ~ 



1 avn + -x'^Vr, ] ba abv + -x^bV 

2 2/ 2 2 



(26) 



Relations (P5|) - (PB| . substituted back into ([22) and ([23]), allow us to exclude the variations of the 
metric and to obtain a system of two coupled equations (each of them of second order) for bF and 
b^. 

By the following substitution, 

c * 

(5F=— ^, (5$ = -, (27) 



we obtain a system of coupled wave equations 

- e"°->^«C + C" + ^"' ^""' C' + ^iC + M^' + M"^ = 0, (28) 

_gao-Xov[, ^ ^n ^ ^2__i!2_,jf' + B^-^ + B^C! + S3C = 0. (29) 

The coefficients A\, A2, A^jBi, B2, and B3 are given in Appendix \K\ If we multiply ([25]) and ([2^ by 
gXo-ao g^jj(j introduce a new radial coordinate 

da;* = — ^ TT^r, (30) 

the system of equations takes the form 

-^ + 4 + M + A.- + A,*.0, (31) 

-ir^ + ^rT + ^i* + ^2T7^ + ^3C = o. (32) 



The coefficients Ai, A2, A3, Si, B2, and B^ are given in Appendix \X\ 

Let us describe the qualitative properties of the wave equations. If we had one wave equation 
(which could be transformed to a stationary Schrodinger-like equation by a proper separation of the 
time and spatial variables) the QNM frequencies would depend strongly on the shape of the potential. 
In our case we have two coupled wave equations and the notion of potential is not so clear. Still, 
the coefficients in front of the zeroth order derivatives of the wave functions, Ai and Bi, respectively, 
determine the properties of the solutions. Thus we introduce two functions [/'■ and C/* which can 
be expressed by Ai and Bi and which we will call potentials of the wave equations for ^ and ^P, 
respectively. As the results show, the presence of unstable modes (solutions divergent with time), 
depends on whether these potentials have a deep enough negative minimum. 

The potentials C/^ and C/* for the solitons are defined as 

U< = _eXo-ao 1^^ __\=_Ia, -^j , [/* = -i?i, (33) 

where C/^ is chosen in such a way that it is regular at the origin (see [llj for a more detailed discussion 
on that definition) and [/* is also finite at a; = 0. The potentials for the black- hole solutions are 
simply 

W^ - -ii, [/* = -Bi, (34) 

where both C/^ and C/* are zero on the event horizon xh- 

It is also important to comment on the asymptotic value of the potentials at infinity. For both 
the black holes and the solitons, U^ tends to zero when a; — J' 00 , but C/* has a nonzero value at 
infinity. The reason is that the scalar field we are considering is massive because of the specific form 
of the potential V{^), defined by eq. ([3]). The mass m of the scalar field can be defined through the 
asymptotic value of the potential C/*, i.e. as 

lim f/* = -m^. (35) 

x—^oo 

Using equation ([23l) it can be easily derived that 

lini (^^] = -m'. (36) 

x^oo \^ 2Na J 

One of the main differences in the time evolution of a massive test scalar field in comparison to 
the massless case, is that the tail is oscillating with the period [32] 

27r , , 

r = — . 37) 

m 

In the limit to — ^ the oscillations of the tail disappear and we are left with the standard power law 
tail. Even though in our problem we are dealing with two wave equations - one for the Skyrme field 
and one for the massive scalar field - it is expected (and confirmed by the numerical results) that the 
tail will again be oscillatory with period ((37|) . 

Another point worth mentioning is the qualitative behavior of the QNM frequencies for the stable 
and the unstable modes. As it is well known, the frequencies of the stable modes are complex, where 
the real part is inversely proportional to the period of the oscillations and the imaginary - to the 
damping time. The picture changes when the modes are unstable. In this case the frequencies are 
purely imaginary, i.e. there is no oscillation and the modes grow exponentially with time |39) , |40) . 

4 Numerical results 

4.1 Solitons 

The background soliton solutions have been obtained in [31 . These solutions are topologically non- 
trivial and the integer n that occurs in the boundary condition for the Skyrme field at the origin 
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Figure 1: The F^^Q{a) and $a;=o(a) phase diagrams for sequences of sohton solutions for n = 1 and 
different values of Dcs ■ 



Fx=o = m: is interpreted as the baryon number. Once n is fixed, the soliton solutions obtained in |31| 
are labeled by the values of the shooting parameters F^^q and $k=Oj where the index {..)x=o refers 
to the value of the function calculated at the origin a; = 0. An example of the F^^Q(a) and '^x=Q{a) 
phase diagrams, presenting sequences of soliton solutions for n = 1, is shown in Fig. [TJ From the 
figure it can be seen that the solutions are divided into two branches - the so-called upper and lower 
branches - and the two branches merge at some critical value of the parameter Ccrit Q Their stability 
is described below. 

We will start with the lower branch of solutions which is stable for ES solitons, i.e. in the case 
without scalar field. The so-called potentials defined by equations (|33)) are given in Fig. [2] for some of 
the soliton solutions which belong to the lower branch in Fig. [TJ As it can be seen the potential C/* is 
positive and cannot lead to instabilities but C/^ is negative near the origin which means that unstable 
modes could exist. The asymptotic value at infinity of C/^ is zero, and for the chosen parameters, C/* 
tends to C/* °° — —0.00346 which means that the mass of the scalar field is to = 0.0589 and the period 
of the oscillation of the tail is T = 107 according to eq. ((37|) . 

We evolve the coupled wave equations (P5)) ~ (|^51) with the appropriate QNM boundary conditions 
- the perturbations should be regular at the origin a; = (i.e. in our case Ca;=o = and ^ x=o — 0) 
and have the form of an outgoing wave at infinity. It turns out that all of the studied solutions which 
belong to the lower branch are stable against the considered perturbations. The time evolution of a 
Gaussian initial perturbation is presented in Fig. [S] The wave form consists of quasinormal oscillations 
in early times and an oscillatory tail for late times, where the period of the tail oscillations is the same 
as the period predicted by eq. (j37p within numerical errors. 

Now let us consider the upper branch of soliton solutions which is unstable for the ES solitons. 
The two potentials [/'» and C/* are shown in Fig. [H Again [/* is positive and V^ is negative near 
the origin but on this branch the negative part is deeper than on the lower branch. This may lead 
to instabilities in the wave equation for C, which will also affect the perturbations of the scalar field 
\1/ through the coupling terms. When we evolve equations (l^51) - (P5|) . it turns out that all of the 
studied upper-branch solutions are unstable. The time evolution of a Gaussian initial perturbation is 
presented in Fig. [S] As it can be seen, the perturbations Q and '^ grow exponentially with time, i.e. 
the QNM frequencies are purely imaginary. 

The calculated frequencies for both of the branches are shown in Fig. [6l The error of the obtained 
frequencies is big and can reach up to 20% in some cases. The reason for this is the more complicated 
wave form which is due to the coupling of the wave equations and the presence of mass term in the 
wave equation for the perturbations of the scalar field ^. Also, the background solutions are known 



*There is a discrete infinite series of copies of these branches corresponding to higher excitations with n > |31| but 
they will not be discussed here since they are energetically unstable 1131 1171 . 
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Figure 2: The potentials U'' and t/* of the lower-branch solitons for Doff = 0.00025 and for several 
values of a. 
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Figure 3: The time evolution of the perturbations C and ^ of a lower-branch soliton solution with 
D^s = 0.00025 and a = 0.2. 

only numerically, which is an additional complication. But even if we take into account the error, the 
numerical values of the QNM frequencies differ significantly from the case without a scalar field [TS] 
(the absolute values of the real and the imaginary parts of the QNM frequencies can be several times 
bigger here than in the ES case). 

The qualitative behavior of the frequencies as we vary the parameter a, which is shown in Fig. [6l 
is the one we expected from the case without a scalar field |T5] . When we increase the value of the 
parameter a, the real uir and the imaginary w/ parts of the frequencies of the stable lower branch 
decrease, while the frequencies w/ of the unstable modes of the upper branch increase (as we already 
said ton = for the unstable upper branch). In the limit a — >■ flcrit (i-e. when we approach the value 
of the parameter a where the two branches merge), the loj of the upper and the lower branches tend 
to zero, i.e. they indicate a stability change. 

When we increase the value of the parameter DeS the errors of the obtained QNM frequencies 
increase mainly because the background solutions become more difficult to obtain and the change in 
the frequencies as we vary Dcs is within numerical errors. 



4.2 Black Holes 

Because of the presence of an event horizon, the Skyrme black holes are topologically trivial. The 
shooting parameters for the background black-hole solutions are the values of the Skyrmion and the 
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Figure 4: The potentials U^ and [/* of the upper-branch soUtons for DefF = 0.00025 and for several 
values of a. 
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Figure 5: The time evolution of the perturbations C, and '^ of an upper- branch soliton solution with 
D^s = 0.00025 and a = 0.2. 



scalar fields at the horizon xh ~ Fh and ^h, respectively. The Fh{xh) and ^h{xh) phase diagrams 
for sequences of black- hole solutions, obtained in [31], are shown in Fig. [T] As we can see, again two 
branches of solutions exist (upper and lower) which merge at some critical value of the radius of the 
horizon rncrit Q The stability of the two branches is described below. 

We will start with the upper branch which is stable for the ES black holes, i.e. in the case without 
a scalar field. The potentials C/^ and t/* for some of the upper-branch solutions are given in Fig. [71 
are given in Fig. [51 Similar to the soliton case, the potential C/* is positive and C/^ has a negative 
minimum near the horizon. We evolve equation (|3i p - (|32p with the standard boundary conditions - 
purely ingoing waves at the horizon and purely outgoing waves at infinity. It turns out that all of the 
studied black holes of the upper branch are stable against the considered perturbations and the time 
evolution of a Gaussian perturbation is shown in Fig. [9l Again, two stages of the time evolution are 
observed - the quasinormal ringing and the oscillatory tail with a period given by eq. p7p . 

The potentials C/^ and t/* for some of the solutions which belong to the lower branch are shown 
in Fig. [TOllj. U is again positive but U'^ has a negative minimum near the horizon which is 



®The upper and the lower branches are defined using the F[j{xfj) diagram. This is obviously different from the soliton 
case, and actually, the upper branch for black holes will have properties (such as stability and finiteness/divergency of 
some of the functions as a — )■ 0) similar to the lower branch for solitons and vice versa. 

^The lower branch is unstable for the ES black holes, i.e. in the case without a scalar field 
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Figure 6: The real (left panel) and the imaginary (right panel) parts of the frequencies as a function 
of the parameter a for the lower and the upper branches of soliton solutions (-Doff — 0.00025). The 
frequencies are obtained using the time evolution. 
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Figure 7: The Fh{xh) and ^h{xh) phase diagrams for sequences of black- hole solutions for D^s 
0.005 and for different values of a. 



much deeper here than for the corresponding black holes of the upper branch, which could lead to 
instabilities. Indeed the time evolution shows that all of the solutions of the lower branch are unstable. 
The logarithm of the wave functions are shown in Fig. 1111 where both C and '^ grow exponentially 
with time. 

The corresponding QNM frequencies for the upper and the lower branches are shown in Fig. 1121 
The error is again big (can reach up to 20%) but we can comment on the qualitative behavior of the 
frequencies. The real and the imaginary parts of the stable upper-branch frequencies decrease when 
we approach the critical value XHcrit where the two branches merge. On the lower branch, which 
is unstable, the imaginary part of the frequencies increases when we increase xh (j-^r = for this 
branch). So, on both branches, in the limit xh — > a^Hcrit, the imaginary parts of the frequencies uJi 
become zero, i.e. a change of stability is observed. 

As discussed in ^31, . the properties of the black- hole solutions when we vary a for fixed xh are 
similar to those of the solitons - two branches of black holes exit which merge at some critical value 
of the parameter Odit. The behavior of the QNM frequencies is also similar to the one shown in Fig. 
[6] for the soliton case. With the increase of the parameter a the real and the imaginary parts of the 
frequencies of the stable branch decreases, and the imaginary part of the frequencies of the unstable 
branch increases. The imaginary parts of both the stable and the unstable branches of black holes 
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Figure 8: The potentials V^ and C/* of the upper-branch black holes for D^s = 0.00025 and for 
a = 0.1. 
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Figure 9: The time evolution of the perturbations C and ^ of an upper-branch black- holes solution 
with Doff = 0.00025, a = 0.1, and r/f = 0.1. 

tend zero when Ocrit is approached. Similar to the soliton case, the change in the frequencies when 
I?cff varies is within numerical errors. 

5 Summary of the results 

As the results indicate, the dilaton does not change the stability of the solutions. Again, the so-called 
lower branch of solitons is stable against radial perturbations while the upper is unstable. For the 
black holes the upper branch is stable and the lower is unstable. The thermodynamical stability 
analysis of the black holes presented in [3T], which is based on the turning point method, is also in 
agreement with the results from the linear stability analysis. 

The modes of the unstable solutions are purely imaginary, i.e ujji — and oj/ < 0. The modes of 
the stable solutions are damped oscillations with uib. 7^ and uji > 0. Hence, at the point of stability 
change both the imaginary part and the real part of the QNM frequencies become zero. The time 
evolution of the black-hole solutions in the vicinity of the point where the two branches merge could 
not be studied accurately, but still the results presented in Figs. [5] and [T^ for the solitons and the 
black holes, respectively, show the expected qualitative behavior. Also, the numerical values of the 
QNM frequencies can be significantly different from the case without a scalar field, i.e. the presence 
of a scalar field significantly alters the spectrum quantitatively. 
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Figure 10: The potentials U'' and [/* of the lower-branch black holes for Doff ~ 0.00025 and for 
a = 0.1. 




ln(v) 




Figure 11: The time evolution of the perturbations (^ and ^ of a lower-branch black-hole solution with 
L'cff = 0.00025, a = 0.1, and r/f = 0.1. 

One more interesting observation can be made. If we consider the evolution of a test scalar field 
on both black-hole and soliton backgrounds, it turns out that all of the modes are damped, i.e. they 
are stable for both branches - the stable and the unstable ones. Thus the time evolution of the test 
scalar field cannot be used to study the stability of the branches. 
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Figure 12: The real (left panel) and the imaginary (right panel) parts of the frequencies as a function 
of the radius of the horizon xh for the lower and the upper branches of black-hole solutions (I?cff = 
0.00025 and a — 0.1). The frequencies are obtained using the time evolution. 



A Coefficients in the wave equations 

The coefficients in the wave equations Eqs. ([25]) and 1^^ are 



A, 



Xo' -ap' uppFf) _ a^^^ 
2 uq X 



2x2' 



-e"°«oFFo' - ^e"°Mo«oFo'' 



2 2 Wo 2 Uq Uq 



,'2 



t'2 



1 Uo2 _ 1 Fq uqfp ^ 1 ^^^ vqff Imq a 

4 Ug 2 Uq 2 Uq 2 Uq x"^ 

^ ^ UQq.FQ 



e'^'uoF, 



'2 
I 



XWUq 



(38) 
(39) 






2 x^/uq ' 2 x^/uq x^^/uq 



1 "OF* j;,'2 1 



2 Xy/^M) 



F^- - -aNe^°^% + aJVe"» ^^^° '^° 



waF_ 
-^e"°i>o*V^i^o + \aNxe'^°^QVQFQ% + J^Fq" 



(40) 



B, 



1 e"Uo*$ Xo'-^o' 1'"0*$j:,'2 



2x 2 iVa;2 



2 iVx2 
xe"yo**o + i 



a; 



So 



NxJu^ 



■Fq 



(41) 
(42) 
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^' = 2 — ^ — ^"-y^^^o'^o + ^^^o'^o 

2N Xyj^ ° 2N x^ ° 2N x^/u^ 

+ - .e VoV^oF, % - - -^$0 + ^ ^^n (43) 



and the coefficients in Eqs. ((3T|) and (|32]) are 



i^ ^ e^o-^Mi, is = e(^°-"°)/2^2, A3 = eX"-"M3, (44) 

ii=eX°""Mi, B2^e^^°-°'°^'^B2, B3 = e^°-"°B3. (45) 
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